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The thermodynamical stability of A1(BH4)3 has been investigated using first-principles 
calculations based on density functional theory. The heats of formation are obtained to be —132 and 
— 131 kj/mol without the zero-point energy corrections for a- and /3-Al(BH~4)3, respectively, which 
are made up of discrete molecular A1(BH4)3 units. It is predicted correctly that the a phase is more 
stable than the phase. The energy difference between the solid phases and the isolated molecule 
is only about 10 kj/mol. An analysis of the electronic structure also suggests the weak interaction 
between A1(BH4)3 molecules in the solid phases. It is confirmed that A1(BH4)3 obeys the linear 
relationship between the heat of formation and the Pauling electronegativity of the cation, which 
has been proposed in our previous study [Nakamori et al, Phys. Rev. B 74, 045126 (2006)]. 

PACS numbers: 65.40.-b, 71.20.Ps, 61.66.Fn 



I. INTRODUCTION 

Metal borohydrides are potential candidates for 
hydrogen storage materials because of their high 
gravimetric density of hydrogen. Among borohydrides, 
alkali borohydrides such as LiBH4 and NaBH4 are 
well known and their properties have been investigated 
moderatelyi ^ViVi 6 ! 7 Since these compounds are 
thermodynamically too stable and desorb hydrogen 
only at elevated temperatures, it is required to decrease 
their hydrogen desorption temperatures for practical 
applications. For LiBELi, the reduction of the enthalpy 
change for the hydrogen desorption reaction has been 
attained by mixing with additives j£*2*i2i whose desorption 
temperatures are lower than that of pure LiBH4 by a 
few hundreds kelvins. 

Although borohydrides composed of non-alkali metals 
can be found in literature, 11 a little is known for 
their properties. In this context, we have recently 
investigated the thermodynamical stability of several 
metal borohydrides, M(BH 4 )„ (M = Li, Na, K, Cu, 
Mg, Zn, Sc, Zr, and Hf), both theoretically and 
experimentally^ It has been found that the stability of 
metal borohydrides shows a good correlation with the 
electronegativity of cations M. We have also determined 
the structural parameters experimentally for Ca(BH4)2 
and its fundamental properties have been predicted 
theoretically>i^ 

Aluminum borohydirde A1(BH 4 )3 is liquid at ambient 
temperatures with the melting point of 209 K. The 
structures of the solid phase have been investigated 
by x-ray diffraction measurements^ Cooling liquid 
A1(BH4)3, the orthorhombic j3 phase was initially grown 
and then the transition to the monoclinic a phase 
occurred at a temperatures in the range 180-195 K. In 
this study, we predict the thermodynamical stability of 
A1(BH4)3 and reexamine the correlation between the 



stability of borohydrides and the cation electronegativity. 

II. METHOD 

The present calculations have been performed 
using the ultrasoft pseudopotential method 1 - based 
on density functional theory. 16 The generalized 
gradient approximation 17 (GGA) is adopted for the 
exchange-correlation energy. The cutoff energies used in 
this study are 15 and 120 hartrees for the pseudowave 
functions and the charge density, respectively. The 
fc-point grids for the Brillouin zone integration are 
generated so as to make the edge lengths of the grid 
elements as close to the target value of 0.08 bohr -1 
as possible. These computational conditions give good 
convergence for the total energy within 1 meV/atom. 
The computational details can be found in Ref. and 
the references therein. 



III. RESULTS AND DISCUSSION 

A. Thermodynamical stability of Al(Bfi4)3 

The heat of formation for the following reaction is 
considered, that is, the formation of A1(BH4)3 from the 
elements: 

A1 + 3B + 6H 2 -> A1(BH 4 ) 3 . (1) 

The heat of formation is estimated from the difference of 
the total energies between the left- and right-hand sides 
of Eq. Q . In this section, we ignore the zero-point energy 
(ZPE) contribution to the heat of formation. 

The crystal structures of a- and /3-Al(BH 4 ) 3 are shown 
in Figs. ( a ) an d (b), respectively. Both phases are 
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(a) a-Al(BH 4 ) 3 




(b) /3-Al(BH 4 ) 3 



FIG. 1: Crystal Structure of (a) a- and (b) /3-Al(BH 4 ) 3 . 
White, gray, and black spheres indicate Al, B, and H atoms, 
respectively. 



made up of discrete molecular A1(BH 4 ) 3 units in which 
three boron atoms form a planar AIB3 skeleton and are 
bonded to Al atom by two bridging hydrogen atoms. 
The packing of the molecules is slightly denser in the a 
phase than in the (3 phase. The structural optimization is 
performed for both phases starting from the experimental 
configurations ^i&isi where the atomic positions as well as 
the lattice vectors are fully relaxed. 

The optimized structural parameters are given in 
Tables HI and ITU for a- and /3-Al(BH 4 ) 3 , respectively. 
The lattice constants are overestimated: The maximum 
deviation from the experimental values is 4 % for the 
a axis in the a phase and 6 % for the b axis in the 
(3 phase, which are larger than the errors of typical 
GGA calculations. The heats of formation are obtained 
as AH = -132 kJ/mol for a-Al(BH 4 ) 3 and AH = 
— 131 kJ/mol for /3-Al(BH 4 ) 3 . It is predicted correctly 
that the a phase is more stable than the (3 phase. 

The relatively large deviations found for the lattice 
constants are probably due to the weak interaction 
between molecular A1(BH 4 ) 3 units in the solid phases. 
The A1(BH 4 ) 3 molecules are expected to be neutral and 
not bonded strongly to each other. In this case, van 
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FIG. 2: Electronic density of states for a-Al(BH 4 ) 3 . The 
origin of energy is set to be the top of valence states. 
The energy positions of the occupied states for molecular 
A1(BH 4 )3 are indicated by solid circles for comparison 
purpose. 



der Walls interactions play an important role, which 
can not be described properly by the current GGA 
functional. To check this problem, we repeat the 
structural optimization, where the lattice constants are 
kept fixed at the experimental values and the atomic 
positions are only relaxed. The obtained heats of 
formation are —128 and —127 kJ/mol for the a and 
(3 phases, respectively. It can be confirmed that the 
relaxations of the lattice vectors give only a minor effect 
for the energetics of A1(BH 4 ) 3 . 

For comparison, the calculation is performed on 
molecular A1(BH 4 ) 3 which has Z? 3 symmetry. A 
hexagonal supercell with lattice constants of a = c — 
13.2 A is used with single T-point sampling. The heat 
of formation is predicted to be AH — — 122 kJ/mol. As 
expected, AH of molecular A1(BH 4 ) 3 is slightly higher 
than those of the solid states by about 10 kJ/mol. This 
small energy difference must be closely related to low 
melting point of A1(BH 4 ) 3 . In Talbe IIIII we compare 
the bond lengths and angles obtained for molecular 
A1(BH 4 ) 3 with those for two solid phases. The geometries 
of A1(BH 4 ) 3 molecules are essentially unchanged when 
forming the solid phases. The agreement with the 
experimental data is good except for some of hydrogen 
related parameters. This is probably caused by the 
experimental difficulty in identifying H positions due to 
their weak x-ray scattering power. 

Figure |3 depicts the electronic density of states for 
a-Al(BH 4 ) 3 . The electronic structure is nonmetallic with 
a calculated gap of 6.0 eV. The occupied states consist of 
several sharp peaks whose energy positions correspond 
well to those of molecular A1(BH 4 ) 3 . The similar 
correspondence can be found for /3-Al(BH 4 ) 3 . These 
also suggest the weak interaction between A1(BH 4 ) 3 
molecules in the solid phases. 
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TABLE I: Optimized structural parameters for q-A1(BH 4 )3. Space group: C2/c (No. 15). a = 22.834, b = 6.176 and 
c = 22.423 A, and (3 = 111.67°. The unit cell contains two independent A1(BH4)3 molecules. The experimental lattice 
parameters are a = 21.917, b = 5.986 and c = 21.787 A, and /3 = 111.90° 



Wyckoff coordinates Wyckoff coordinates 





position 


X 


y 


z 




position 


X 


y 


z 


Al 


8f 


0.3797 


0.5943 


0.8366 


Al 


8f 


0.3659 


0.1127 


0.5947 


Bl 


8f 


0.3205 


0.3121 


0.8239 


Bl 


8f 


0.3626 


-0.1633 


0.5349 


Hl„ 


8f 


0.3809 


0.3078 


0.8439 


Hl a 


8f 


0.3860 


-0.1645 


0.5965 


Hl 6 


8f 


0.3002 


0.5071 


0.8128 


HI,, 


8f 


0.3437 


0.0276 


0.5138 


Hl c 


8f 


0.3067 


0.2479 


0.8677 


Hl c 


8f 


0.3163 


-0.2759 


0.5191 


Hl d 


8f 


0.3019 


0.2214 


0.7725 


Hid 


8f 


0.4058 


-0.2032 


0.5180 


B2 


8f 


0.3899 


0.7542 


0.7555 


B2 


8f 


0.2793 


0.2182 


0.6028 


H2 a 


8f 


0.4165 


0.5751 


0.7794 


H2 a 


8f 


0.3147 


0.0676 


0.6368 


H2 6 


8f 


0.3516 


0.8141 


0.7820 


H2 6 


8f 


0.3067 


0.3176 


0.5703 


H2 C 


8f 


0.3575 


0.7139 


0.7005 


H2 C 


8f 


0.2749 


0.3494 


0.6407 


H2 d 


8f 


0.4309 


0.8894 


0.7685 


H2d 


8f 


0.2332 


0.1306 


0.5652 


B3 


8f 


0.4298 


0.7324 


0.9297 


B3 


8f 


0.4542 


0.2773 


0.6427 


H3 a 


8f 


0.3703 


0.7520 


0.8985 


H3 a 


8f 


0.4187 


0.3056 


0.5838 


H3 6 


8f 


0.4541 


0.6038 


0.9002 


H3i, 


8f 


0.4282 


0.1495 


0.6705 


H3 C 


8f 


0.4349 


0.6401 


0.9787 


H3 C 


8f 


0.4573 


0.4508 


0.6683 


H3d 


8f 


0.4531 


0.9080 


0.9304 


H3d 


8f 


0.5010 


0.1867 


0.6426 



TABLE II: Optimized structural parameters for /3-Al(BH4)3. 
Space group: Pnalx (No. 33). a = 18.649, b = 6.488 and 
c = 6.389 A. The unit cell contains one independent A1(BH4)3 
molecule. The experimental lattice constants are a — 18.021, 
b = 6.138 and c = 6.199 Ap±^±2 



Wyckoff coordinates 





position 
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y 
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Al 


4a 


0.8703 


0.1558 


0.2098 


Bl 


4a 


0.7800 


0.0057 


0.0633 


Hl„ 


4a 


0.8456 


-0.0341 


0.0213 


HU 


4a 


0.7751 


0.1384 


0.2112 


Hlc 


4a 


0.7552 


-0.1515 


0.1331 


Hid 


4a 


0.7551 


0.0858 


-0.0906 


B2 


4a 


0.9168 


0.0183 


0.4855 


H2 a 


4a 


0.9353 


-0.0248 


0.2979 


H2 6 


4a 


0.8712 


0.1661 


0.4885 


H2 C 


4a 


0.9700 


0.0865 


0.5674 


H2d 


4a 


0.8870 


-0.1319 


0.5547 


B3 


4a 


0.9115 


0.4349 


0.0722 


H3 a 


4a 


0.8623 


0.4281 


0.2121 


H3 6 


4a 


0.9337 


0.2529 


0.0277 


H3 C 


4a 


0.9619 


0.5214 


0.1520 


H3d 


4a 


0.8837 


0.5009 


-0.0844 



B. Correlation between stability of borohydrides 
and cation electronegativity 

In Ref. IT3. we have found a good correlation 
between the heat of formation AH and the Pauling 
electronegativity of the cation xp> when AH is 
normalized by the number of BH4 complexes in the 
formula unit. This correlation can be represented by the 
linear relationship, AH — 248. 7\p — 390.8 in the unit of 
kJ/mol BH4. In this section, we refine this relation using 
the additional results for a-Al(BH4)3 and Ca(BH4)2>^ 

In order to estimate the ZPE contribution to AH 
for Al(Bfi4)3, the normal eigenmode frequencies are 
calculated for molecular Al(Bfi4)3. The weak interaction 
between A1(BH4)3 molecules in the solid phases justifies 
this treatment. The result is shown in Fig. |3 where 
the obtained frequencies agree fairly well with the 
results of the experiments and the quantum chemistry 
calculation^ The ZPE contribution to AH is predicted 
to be 36 kJ/mol BH4. This value is quite close to 
the approximated value of 33 kJ/mol BH4 adopted in 
Ref. Il2l which has been obtained using the molecular 
approximation for a [BH4]~ anion. To keep the 
consistency with the previous treatment, we decide to 
use the approximated value of 33 kJ/mol BH4, where the 
normalized heats of formation with the ZPE correction 
become —11 and —155 kJ/mol BH4 for A1(BH4)3 and 
Ca(BH4) 2 , respectively. 

In Fig. 21 the normalized heat of formation AH 
with the ZPE correction as a function of the Pauling 
electronegativity of the cation xp is plotted^ It can 
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TABLE III: Comparison of the bond lengths d (A) and angles 9 (°) obtained for molecular A1(BH4)3 and two solid phases. 
For d(Al-H) and #(H-A1-H), H atoms forming the bridge bond are only cosidered. 







Calc. 




Expt.° 




molecule 


a phase 


(3 phase 


a phase 


(3 phase 


d(Al-B) 


2.16 


2.15-2.16 


2.15-2.16 


2.10-2.14 


2.10-2.13 


d(Al-H) 


1.77 


1.77-1.78 


1.77-1.78 


1.73-1.76 


1.68-1.75 


d(B-H) 


1.20-1.29 


1.20-1.29 


1.21-1.28 


0.99-1.14 


0.99-1.14 


(9(B-A1-B) 


120 


117-122 


119-122 


119-121 


119-123 


6>(H-A1-H) 


73 


73 


73 


65 


63-65 


(9(H-B-H) 


106-122 


105-121 


105-121 


103-121 


106-130 
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FIG. 3: Normal eigenmode frequencies for molecular 
A1(BH4)3. The calculated frequencies are indicated by 
vertical bars and the Gaussian broadening with a width of 
30 cm -1 is used. 

be found that the results for A1(BH4)3 and Ca(BH4)2 
also obey the linear relationship. The least square fitting 
yields 

AF = 253.6xp - 398.0, (2) 

with an absolute mean error of 9.6 kJ/mol BH4. The 
addition of two new data points causes only a little 
change in the coefficients for the linear relation. The 
present study provides the additional support for the 
linear correlation between AH and xp- The Pauling 
electronegativity of the cation \p 1S a good indicator to 
estimate the thermodynamical stability of borohydrides. 

IV. CONCLUSIONS 

In this study, we have predicted the thermodynamical 
stability of A1(BH4)3 and reexamined the correlation 
between the stability of borohydrides and the cation 
electronegativity. 

The heats of formation are obtained as —132 and 
— 131 kJ/mol without the zero-point energy corrections 
for a- and /3-Al(BH4)3, respectively, which are made up 
of discrete molecular A1(BH4)3 units. It is correctly 




FIG. 4: Relation between the heat of formation AH and 
the Pauling electronegativity of the cation, \p. The straight 
line indicates the result of the least square fitting, AH — 
253. 6xp — 398.0. The zero-point energy contributions to AH 
are approximately taken into consideration (see the text). 

predicted that the a phase is more stable than the (3 
phase. The energy difference between the solid phases 
and the isolated molecule is small and the interaction 
between the molecules in the solid phases is expected to 
be weak. This is most likely related to low melting point 
ofAl(BH 4 ) 3 . 

The linear correlation between the heat of formation 
and the Pauling electronegativity of the cation, which 
has been proposed in the previous study^ is also held 
for A1(BH 4 ) 3 as well as for Ca(BH 4 ) 2 . The Pauling 
electronegativity of the cation is a good indicator to 
estimate the thermodynamical stability of borohydrides. 
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